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Abstract 

By means of direct canonical transformations, the Poincare- 
Von Zeipel perturbation method may be modified into a more useful form<> 
Teo particular variants of this approach are examined and explicit- 
formulas are given for carrying them out to at least the 4-th order. 
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1. INTRODUCTION 

In this work, a modification of the Hamilton-Jacobi equation is developed, 
utilizing the direct formulation of canonical transformations. The resulting 
perturbation scheme is easy to apply and may be related either to the standard 
Poiucare-Von Zeipel method or to Hori's expansion, depending on the choice of 
certain expressions entering it. 

The solution of a mechanical problem by the Hamilton-Jacobi equation seeks 
a canonical transformation to new variables, in which the complexity of the 
problem is reduced. In particular, let us consider a perturbed periodic motion, 
with p 1 the action variable of the unperturbed motion and with a time-inde- 
pendent Hamiltonian 

H = p, + 12 £ k K^(p, q) (l) 

A k*l " ' 

Then by the Poincare-Von Zeipel method (Poincare, 1895; Von Zeipel, 1916; 
Corben and Stehle, I960) one seeks a near-identity transformation to new variables 
(P, £) , given by the generating function 

* H P i \ + 

via 

^<T/'DP i = 

^/'Dq i - P A 

In this method, one uses the fact that the new Hamiltonian equals the 
old one, substituting (5) to obtain 

H(T)r/T>q, a ) . H *(£, W/T>£) U) 


Hi t k <3^(P, q) (2) 

n — 
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When this ia expanded in orders of £ , a recursive scheme is obtained and both 

# 

H and (T may be derived order by order, subject to the condition that the trans- 

* 

formed angle variable is absent from the new Hamiltonian H . When this is 

accomplished, not only has Q 1 been removed from the problem, but its conjugate 

has become a constant of the motion, so that the motion represented by H has 
two fewer variables to contend with. 

Now a near-identity transformation may be represented in other ways than by 
equations (2)-(3) ; these ways may offer more convenience, since C T in (2) depends 
on both old and new variables and therefore the relations (3) require further 
unt angling to be useful. In particular, we may consider a direct canonical trans- 
formation 

* - y + ZZt* $ lk) ly) (5) 

k=l ' 

where we will use the notation 



Then, if one substitutes (3) into the basic relation 
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H(jr) = H {%) (y) 

and expands order by order , the appropriate transformation (5) nay be derived. 

Horl (1966) used this approach, with (9) given through Lie's formulation of 

canonical transformations (Lie, 1888). Here a somewhat different approach (Stern, 

1970) will be used, which can be made equivalent either to Hori's or to the Poincare- 

(k) 

Von Zeipel method, depending on one's choice of the expression f , defined later. 

2. NOTATION 


If the position vector y in phase space is defined as in (6), let its 
conjugate y be defined 


from which, by (7) 


y = ( q* -P ) 



( 10 ) 


Also, if v operates in y-space, we can define there a conjugate gradient 
operator V , the i-th component of which is . 

To express functions of z in terms of £ , we shall vise Taylor expansion 
operators T^ (Musen, 1969)# with operation denoted (for clarity) by an asterisks 

G(z) = G(y + ££ k $ (k) ) 

a exp( S £ k ^ k ^»V) * G(y) 
k=l 

= Z6 k T^* G(y) (12) 

ksO 

(k) 

The T x may be obtained by expanding the exponential and regrouping terms 
according to their order in £ ; the first few of them are 
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T (0) 

® 1 



T (!) 

- 






i (13) 

t (2) 

_ ^(2). y + VV 


1 

t (3) 

*** 

+ (1/6) vw 

and 

so forth. In general one can write 



T (k) 

= $ (k >.V + N^ k ^ 


(14) 


(k) 

where N' ' is an operator involving V at least twice. 

3* THE TRANSFORMATION 

this 

The relation (3) holds for any near-identity transformation* If ^is to 

v(k) 

he a ftftH transformation* $ must satisfy additional requirements* and 

it may he shown (Stern* 1970) that these always have the general form 

$00 . fW(£) + V% (k) (15) 

where ff k \ ^ ) are expressions involving lower orders and is the k-th 

order of a generalization of the generating function. Many alternative choices of 
f^ k ) are possible (one can often construct lower order expressions that have the form 

of a conjugate gradient and -add-thea-to f^ ) hut we shall he concerned with two in 

particular. One may choose 

4 k) - - £ s (,) * o* (lw,) /-5? t > 


(16) 
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where S^ are operators resembling the of (13) but with 

everywhere replaced by 

.TT (j) - (^ j) . ... 0, ... 0 ) 07) 

that is, with vectors in which the momentum-like parts are retained and the 
coordinate-like parts are replaced by zeros. It may then be shown (see Stern, . 1970 . 
where this choice of f^ is denoted by ) that in this case 

* (T (k) ( 18 ) 

with CT W the k-th order "conventional" generating function of equation (2) 
which describes the same transformation. Since contains orders of ^ and 

of X lower than the k-th, equation (16) must be derived recursively order by order, 
Alternatively, f£ k ^ may be derived without breaking up i»to 

momentum-like and coordinate-like parts. The general method for doing this (stern, 
1970) is too lengthy to be described here, and we shall merely give the lowest 
orders of the results 

= 0 

. i 5 (1) - v 1 (1) 

= $ (2) * V5 (l) V(19) 

+ i [(| (1) *v^ (1) )-v| (2) - | (2) *V(| (l) .V^ (l) )j 


f d) 

f (2) 

f (3) 

f U) 
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With this choice it may be shown (Stem, 1970) that 

^< k > . . S W (20) 

Mh.r* S^ k ' is the k-th ord«r generating function defined by Horl (1966). 

4« THE PERTURBATION METHOD 

Expanding ( 9 ), and using (l) and (12;, gives 

y, ♦ Z£ k H U) (y) - E£ k H*W(i) 

A kal k=0 

- Z£ k H* (k) (y ♦2> k 5 (k) ) 
k=0 

- Ee k C iHi'Wffl (21) 

kssO m=0 

This may be equated order by order, the zeroth order being simply 

H* (< %) - H (0) (y) - y x (22) 

In the 0( £ k ) relation we separate the terms with m equal to 0 and k 
from the rest and use (14) and (22). This gives 

H (k) (y) = H* (k) (y) ♦ £ (k) - Vy L . » (k) . yj + 

+ 7Z l (o) . H*( k "“ ) (y) (23) 

Bel 

In the last equation, may be dropped, since it contains V twice or 

more and its aotion on y } thus yields zero. Also, by (ll) 


a 


£^•7 3^ b a + l>% ik) /bY n ( 24 ) 

The final result, which is the basic relation of our perturbation scheme, may 
thus be written 

- A (k) <?> (26) 

where 

f\(*\ y) » H^(y) - fj k \5) - 2U T ^* H*^ k “*\y) (26) 

nal 

If the variable y enters only as an angle *ith basic period unity, 
is clearly periodic, since any nonperiodic part of ~)C is 
removed by the differentiation. If one defines an averaging operator 
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<A W > - 1 A (k) -»„ 

(27) 


0 


then clearly 

H*W») „ 4 A (k) > 

(28) 

x (k) 

. j” ( a (! ^ - <A (k ^> ) ay; + M^Cy) 

(29) 


0 


where M^(y) is an arbitrary function, representing the fact that when y^ and 

v are el imina ted, the other variables may also undergo an arbitrary near-identity 
'n 

transformation among themselves. If no other considerations exist, it may be set 
equal to zero. 
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Suppose the calculation has been carried out up to and including order (k-l). 
We now form A by (26), derive H*^ and by (28) and (29), 

and then obtain from (15), using the appropriate choioe of f^) , The 

derivation is now complete to the k~th order. 
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